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ABSTRACT 


The  pressure  on  an  infinite  sinusoidal  wall  in  steady 
two-dimensional  supersonic  flow  of  an  inviscidgas  in  vibrational 
or  chemical  nonequilibrium  is  obtained  by  the  linearized  method 
of  characteristics  and  compared  to  the  analytical  solution  of  the 
linearized  equations  previously  obtained  by  Vincenti.  The  system 
of  "frozen"  characteristics  is  used.  At  near -equilibrium  con¬ 
ditions  the  error  in  the  method  of  characteristic  s  is  due  entirely 
to  the  error  in  integrating  the  rate  equation  along  the  stream¬ 
lines  and  can  be  made  arbitrarily  small  by  decreasing  the  size 
of  the  characteristic  net.  There  appears  to  be  no  approximate 
criteria  for  reducing  the  error  at  near -equilibrium  conditions 
by  switching  over  to  the  system  of  equilibrium  characteristics 
which  are  applicable  to  the  entire  flow  field. 
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1.  INTRODUCTION 


The  phenomena  of  wave  propagation  and  the  existence  of  characteristic 
surfaces  in  nonequilibrium  flow  have  been  studied  by  several  investigators 
(e.  g.,  References  1  and  2).  It  has  been  well  established  that  character¬ 
istics  exist  for  supersonic  steady  flow  and  are  determined  by  the  frozen 
speed  of  sound  for  all  finite  values  of  the  reaction  rate  regardless  of  how 
large  the  reaction  rate  becomes.  When  the  reaction  rate  becomes  infinite, 
however,  the  characteristics  are  determined  by  the  equilibrium  speed  of 
sound.  This  discontinuous  change  in  the  identity  of  the  characteristics  is 
related  to  a  reduction  in  the  order  of  the  governing  differential  equations 
for  an  infinite  reaction  rate.  It  has  been  shown,  however,  that  the  ana¬ 
lytic  solution  for  a  flow  field  is  a  continuous  function  of  reaction  rate 
(e.  g.,  References  2  and  3).  Thus  a  near  equilibrium  solution  approaches 
the  usual  equilibrium  solution  as  the  reaction  rate  approaches  infinity. 

If  it  is  desired  to  set  up  a  calculation  procedure  based  upon  the 
method  of  characteristics  that  will  be  valid  for  all  reaction  rates  the  sys¬ 
tem  of  frozen  characteristics  must  be  used.  Since  only  finite  reaction 
rates  are  physically  possible  anyway,  such  a  procedure  would  appear  to  be 
quite  satisfactory.  However,  the  characteristic  equations  can  in  general 
be  solved  only  by  use  of  finite  difference  methods,  and  certain  errors  are 
therefore  introduced.  In  the  usual  calculations  for  perfect  or  equilibrium 
gases  these  errors  are  known  to  be  dependent  on  the  size  of  the  character¬ 
istic  net,  Ax/t  ,  where  i  is  a  characteristic  dimension  of  the  flow  field. 
The  introduction  of  a  new  dimension  into  the  problem,  namely  the  relaxa¬ 
tion  length  intuitively  suggests  that  the  error  will  depend  also  on 

the  parameter  Ax/t^U^. 

Practical  considerations  in  calculating  a  flow  field  place  a  certain 
lower  limit  on  Ax/i.  Thus  as  the  reaction  rate  becomes  large  for  a  fixed 
value  of  Ax/i  the  parameter  Ax/x^U^  becomes  large  and  suggests  large 
errors.  For  near  equilibrium  flow  it  seems  likely  that  there  is  a  certain 
value  of  the  reaction  rate  above  which  the  error  could  be  reduced  by 
switching  over  to  the  system  of  equilibrium  characteristics. 
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The  work  described  in  this  report  was  undertaken  in  an  attempt  to 
provide  a  better  understanding  of  the  errors  associated  with  the  calculation 
of  a  non  equilibrium  flow  field  by  the  method  of  characteristics.  In  partic¬ 
ular,  criteria  were  sought  that  will  predict  at  what  large  value  of  reaction 
rate  the  equilibrium  characteristics  should  be  used.  The  analytic  solution 
obtained  by  Vincenti,  Reference  3,  for  the  linearized  non  equilibrium  flow 
past  an  infinite  sinusoidal  wall  was  compared  with  results  obtained  by  the 
method  of  characteristics  solution  of  the  same  linearized  flow  equations. 
Vincenti's  analytic  solution  has  the  advantage  of  being  relatively  simple  so 
that  the  algebra  involved  is  not  prohibitive.  The  linearized  equations  offer 
another  advantage  in  that  for  frozen  flow  the  solution  obtained  by  the  method 
of  characteristics  is  exact  for  all  net  size.  All  errors  can  therefore  be 
attributed  directly  to  the  nonequilibrium  aspects  of  the  flow. 
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2.  LINEARIZED  EQUATIONS  AND  SOLUTION  FOR  A  WAVY  WALL 


The  linearized  equations  of  motion  for  a  gas  subject  to  a  single  re¬ 
laxation  process  were  derived  by  Vincenti  in  Reference  3.  It  was  assumed 
that  the  flow  field  is  described  by  a  perturbation  of  a  uniform  parallel  flow 
with  velocity  U  in  the  x  direction  such  that  u  =  U__  +  u',  v  =  v', 
p  =  p^  +  p',  etc.  It  was  also  assumed  that  the  uniform  flow  is  in  equilib¬ 
rium.  The  resulting  equations  are: 


Continuity: 


■^00 


feu'  .  ev'v 


+  U  4^1  = 
00  8x 


(1) 


X  Momentum: 


■^oo 


U  ^ 
oo  8x 


0 


(2) 


y  Momentum: 


■^00 


U 

00 


8v'  ^ 

8x  8y 


0 


(3) 


Energy: 


u  4^:  + 

00  8x 


(4) 


Equation  of  State: 


dh'  = 


dq' 


(5) 


Rate  Equation: 


8x 


X 


oo  00 


Equilibrium  Relation: 


(6) 


(7) 


The  coordinate  system  used  is  shown  in  Figure  1. 
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The  quantity  q  describes  the  vibrational  or  chemical  state  of  the  gas. 
For  vibrational  relaxation  q  is  the  vibrational  temperature,  and  for  chem¬ 
ical  relaxation  q  is  the  dissociation  mass  fraction.  The  quantity  q  is  a 
fictituous  local  equilibrium  value  of  q  determined  by  the  nonequilibrium 
values  of  p  and  p. 

Vincenti  analyzed  the  flow  past  an  infinite  sinusoidal  wall  by  intro¬ 
ducing  a  velocity  potential  and  solving  the  resulting  third  order  linear 
potential  equation.  Details  of  the  method  and  results  are  contained  in 
Reference  3.  The  perturbation  velocities  are: 


<  (6  .in  2,.  CO.  2. 

00  6  +  V  '  ' 


(8) 


ir  -  ^"7 

00 


e  cos  ZiT  ^ 


(9) 


The  parameters  6  and  X  are  given  by 


Disturbances  in  the  flow  field  are  propagated  along  the  lines  (dy/dx)  =  l/X, 
and  5  is  a  measure  of  the  rate  of  decay  of  these  disturbances.  In  the  equi¬ 
librium  and  frozen  limit  for  supersonic  flow  these  parameters  become: 

Equilibrium  Frozen 

6 
\ 
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3.  SOLUTION  BY  METHOD  OF  CHARACTERISTICS 


The  characteristic  relations  for  the  linearized  Equations  (1)  through 
(7)  were  obtained  by  the  usual  method  of  solving  for  a  particular  derivative, 
such  as  8u'/8x  and  setting  the  determinants  in  the  numerator  and  denomi¬ 
nator  equal  to  zero.  For  all  finite  reaction  rates,  0  ^  i  /t  U  oo,  the 
characteristics  are: 

along 

(13) 


“F 


-  o' 


*  PooHx>  ®F  ^  7^ 

00  00 


^dx=  0 


along 


(14) 


where 


(15) 


(16) 


Thus  the  characteristic  direction  is  determined  by  the  frozen  speed  of 
sound  for  all  finite  reaction  rates.  In  the  linearized  approximation  the 
lines  (dy/dx)  =  0  are  streamlines  and  Equation  (14)  is  just  the  rate  equa¬ 
tion  following  a  fluid  particle. 
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For  flow  that  is  everywhere  identically  in  equilibrium,  ^/■’’qq^oo  = 
the  characteristic  relations  can  be  obtained  by  substituting  q  =  q'  into 
Equations  (1)  through  (7)  and  proceeding  as  before.  The  results  are: 

along 

^  =  i  tan  a 
dx  e 

dp'  ±  p  U  tan  a  dv'  =  0  (17) 

oo  oo  e 


where 


(18) 


The  discontinuous  change  from  the  frozen  to  the  equilibrium  charac¬ 
teristics  as  ■^/'^QQ^QQ  — ^  ^  brings  up  the  question  of  what  errors  are  intro¬ 
duced  in  calculating  near  equilibrium  flow  using  the  frozen  characteristics. 
In  equilibrium  gas  calculations  the  error  which  results  from  solving  the 
characteristic  equations  by  finite  difference  methods  is  known  to  depend  on 
the  size  of  the  characteristic  net,  Ax/i.  The  introduction  of  a  new  charac¬ 
teristic  length  into  the  problem  (i.  e.,  the  relaxation  length,  t^U^) 
suggests  that  the  error  will  now  depend  also  on  the  parameter  Ax/x^U^. 

To  study  this  error,  the  pressure  at  a  wall  point  as  calculated  by  the 

method  of  characteristics,  p^,  is  compared  to  Vincenti's  analytic  solution, 

p  .  A  convenient  parameter  for  expressing  the  error  is  defined  as: 

SI 

AJL  =  _  -  (19) 

Ap  (*^e  PFjmax  along  wall 

where  p^  is  the  pressure  at  the  wall  point  3  (see  Figure  1).  This  param¬ 
eter  is  calculated  as  a  function  of  net  size,  Ax/i  ,  and  reaction  rate, 
i  /t  U  .  All  quantities  along  the  vertical  line  x  =  ^  and  the  reference 

♦  >!c 

pressure  Ap  are  taken  from  the  analytic  solution,  and  Ap  is  given  by: 


Ap 


=  Ztt  p  U 
'oo  oo 


7<.“ 


-  tan 


“e) 


(20) 
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Applying  Equation  (13)  in  finite  difference  from  along  the  character¬ 
istic  1-3  gives: 

P3  =  Pi'  +  PqoHx)  ‘^F  ^^3  '  '^l')  r  B^l'  ■  'll')  +  <^3  "  *133 

00  00 

The  term  q^'  can  be  eliminated  by  noting  that: 

qj’  =  ni'  +  (^13'  -  qi')  (22) 


The  differential  form  of  Equations  (2)  and  (4)  can  be  written  as: 

p  U  du'  +  dp'  =  0 
'oo  00  ^ 

U  du'  +  dh'  =  0 
00 

Equations  (23)  and  (24)  can  be  combined  with  (5)  and  (7)  to  give: 


d5  '  =  Ddp '  +  Edq ' 


where 


D  = 


15  13  iii 

IS  +_8JEL  .  i£^ 

8p  ^ 

Pap  ap  j 


00 


E  = 


8q  8h 
Tp  8q 
8h 

L  ^  J 


00 


Equation. (22)  can  now  be  written  as: 


**3  "  'll'  °*P3  •  Pi')  ^(*13  "  ^1') 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


To  obtain  the  pressure  at  3  from  Equations  (21)  and  (28),  the  rate 
equation  [Equation  (14^  must  be  integrated  along  the  wall  streamline  from 
2  to  3  to  obtain  q^'.  It  is  convenient  and  instructive  to  assume  that  this 
integration  introduces  an  error  Aq'  in  q^'  defined  by: 

^3  =  (<13  )a  +  (29) 
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Substituting  Equations  (28)  and  (29)  into  (21)  gives  the  error  in  pressure  as: 


Ap 


[2^^'  -q^'  -Dp^- 

-Eq^'+(E-  1)| 

Ax 

T  U 

00  00 

\  A  *  1 

Ap 

A®  Ax  ^ 

i  T  u  1 

^  oo  ooj 

Ap 


00  oo  Aq’ 

“ — I55~a5E —  7^ 

‘  -  T-  nr-  ^0 

GO  OO 


(30) 


where  v^'  can  be  evaluated  using  the  linearized  boundary  condition  at  the  wall: 


(V)'  =  2ir  U  .  -  cos 


'wall 


00 


■>  X 

2nj 


(31) 


Thus  the  error  consists  of  two  distinct  parts.  The  term  in  ^  ^  in 
Equation  (30)  is  the  error  which  results  from  the  averaging  of  quantities 
in  the  finite  difference  relation  along  the  inclined  characteristic  1-3.  The 
remaining  term  in  Equation  (30)  is  the  error  in  pressure  associated  with 
the  error  in  integrating  the  rate  equation  along  the  streamline  2-3.  We 
denote  these  in  errors  I  and  II  respectively  so  that: 


(32) 


Vincenti's  analytic  solution  is  used  to  evaluate  the  term  (Ap/Ap  )j . 
The  resulting  complicated  expression  can  be  simplified  by  expanding  it  in 
a  Taylor  series  about  Ax/i  =  0  giving: 


where 


K 


AD  i 
4^  '  ~r  THJ 


GO  00 


(33) 


(34) 
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The  parameters  and  were  determined  for  two  values  of  4>: 
at  <|>  =  0: 


o  3  ^  2  - 

Ztt  tan  Up  o 


X 


K.  =  »i.  . - r-  \  3(\  tan  a  —  )  +  4\  tan  a-,  + - - 

1  3(tan  ap  -  tan  F  F  t^n  ap)‘ 


(35) 


^2  “  3(tan 


*  ^  r- 

TT  tan  a^,  j  a  ^  2 

^  r  <  (\  tan  Qp)  +  (X  tan  Op)  -  2(X  tan  Up) 
F  V 


2X  tan  Op  +  1  +  (X  tan  Up) 


'} 


at  4>  =  f  /4: 

K.  = 


2it^  J  .  .4  ,,,  .  .3  X^ 

1  =  yxftan  a-  t^-a;T  V  °F>  ’  ‘an  Up) 


(36) 


+  ^-5 - y  +  2X  tan  Q  ^ 

6^  +  X^ 


(37) 


K 


2  =  — 2 - - -  /4(X  tan  Op)^  +  3(X  tan  Up)"^  j 

3X^(tan  -  tan  a_^)  1  6  +  X 


-  4(X  tan  ap)^  -  2(X  tan  Up)^  -  (38) 

6  +  X 

In  both  the  equilibrium  and  frozen  limit  the  quantity  (Ap/Ap* )j  goes  to  zero. 

The  error  in  calculating  the  pressure  by  the  linearized  method  of  charac¬ 
teristics  for  near  equilibrium  flow  is  therefore  due  entirely  to  the  error 
in  integrating  the  rate  equation  along  the  streamline. 

In  evaluating  the  term  (Ap/Ap  )j^  it  is  convenient  to  introduce  the 
pressure  coefficient. 


P  = 


P  -  R 


oo 


1/2  p  U  ' 

'  ^X>  00 


such  that: 


where: 


K 


Ax 

9  i 


Aq' 


‘  +  K4  I’^o. 


Ko  =  A(E  -  1) 


P  u 
•oo  00 


00  I 


00  CD 


(39) 


(40) 


A  first  approximation  for  q^'  is: 


Ax 


(41) 


The  error  in  q^' 
Ax  =  0: 


3  -t.  V“*/2 

can  be  found  by  expanding  q'(x)  in  a  Taylor  series  about 


Ax  + 


(42) 


Subtracting  Equation  (42)  from  (41)  gives: 

..(i) 


Aq' 


■  ■  i&), 


A  2  d 

Ax 


-  IT  (^) 


dx 


dx 


(43) 


The  derivatives  in  Equation  (43)  are  evaluated  using  Vincenti's  solution  at 
X  =  ^.  This  results  in: 


(44) 


T  AP 
00 
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where: 


at  4>  =  0: 


The  error  in  q^'  resulting  from  integration  of  the  rate  equation  along 
the  streamline  can  be  reduced  by  use  of  an  iterative  procedure  as  follows: 
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(52) 


(  ) 


and  expand  dq'/dx  in  a  Taylor  series  about  Ax  =  O': 

(»).  -fe-ft) 


Ax^  + 


Solving  Equation  (52)  for  the  quantity  [[  3  substituting  this  into 
Equation  (51)  gives: 


(‘IsOa  =  *12  +  I  ^  i  (^)3  ■  ■re  (^) 


Ax3  + 


(53) 


Subtracting  Equation  (53)  from  (50)  gives: 


'dq'V”'**  ^dq'1 


Ax^  + 


(54) 


The  rate  equation.  Equation  (14),  can  be  used  to  express  Equation  (54). in 
terms  of  the  errors  in  q  '  and  q*  resulting  in: 


Ax^  + 


(55) 


The  error  in  q'  after  (n  -  1)  iterations  is  obtained  from  Equation  (25): 


^i(n-l)  _  +  EAq'^“"^^ 


Aq 


EAq' 


(56) 


Substituting  Equation  (56)  into  (55)  and  using  the  parameters  introduced  in 
Equation  (44)  for  ^d^q'/dx^j^  ^d^q'/dx^^2  gives: 


,(n) 


00 


where: 


./.  \(n-i)  /A  \  A  i(>^“l)  /A  \3  2  /a  \4 

'  00^'^  (57) 


Dp  U 

„  _  .  'oo  00 
*^5  ■  “if 


i 

T.ir 


(58) 


00  00  00 


1 
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(59) 


! 

i 


(  I 


“6  =  I  <>=  -  ‘I  tV 

00  00 


Equation  (57)  can  now  be  used  with  Equations  (39)  and  (44)  to  find 
(Ap/Ap*)'^^  after  an  arbitrary  number  of  iterations.  The  total  error  in 
pressure  after  each  iteration,  (ap/Ap*)^’^^  is  given  by: 


(60) 


(  I 


,1 
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4.  NUMERICAL  RESULTS  AND  CONCLUSIONS 


(  I 


( 


1 


‘  ?  I 

i 

i 


The  results  of  Section  3  are  applicable  to  a  gas  in  which  either 
vibrational  or  chemical  relaxation  occurs.  Numerical  values  of  the  error 
parameter  were  calculated  for  a  diatomic  gas  subject  to  vibrational  relax¬ 
ation  for  which  q  is  taken  to  be  the  vibrational  temperature,  T^^.  The 
equation  of  state  and  the  enthalpy  are  given  by; 

p=pRT  (61) 

h  =  c  T  + 

P 

where  c^  is  the  vibrational  heat  capacity.  The  parameters  in  Section  3 
become: 


/ 


c.  d  T. 


(62) 


IT  1^2 

K4  =  5tan  Uj,  - 


*^ioo 


M 


00  00  00 


tr  _  1  \  #  2  i 

^5 '  VS  Tvr 

00  00  oo 


(63) 

(64) 


^6 


(65) 


where  c^^  is  Cj(Tj)  evaluated  at  the  free-stream  equilibrium  temperature, 

T  .  Results  are  shown  in  Figures  2  and  3  with  initial  conditions  taken 
oo 

along  ^  =  0  and  ^  =  i/4  for  various  numbers  of  iterations.  The  limit 
error  curves  obtained  after  a  large  number  of  iterations  are  replotted 
in  Figure  4  for  comparison. 

The  calculation  procedure  using  the  frozen  characteristics  is  seen 
to  be  convergent  even  in  the  limit  of  equilibrium  flow.  The  error  con¬ 
verges  to  a  limit  error  curve  after  a  number  of  iterations  and  approaches 
zero  with  decreasing  net  size.  Since  all  errors  are  of  order  (Ax/i  )^  or 
higher,  the  error  in  calculating  a  flow  field  for  all  values  of  f  and 

n  can  be  made  arbitrarily  small  by  decreasing  the  size  of  the  character¬ 
istic  net. 
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It  is  noted  that  for  i/x  U  — ♦  oo  only  a  single  iteration  is  required 

00  00 

to  reach  the  limit  error  curve.  This  can  be  explained  by  recalling  that  in 
the  equilibrium  limit  no  error  is  introduced  in  the  calculation  along  the 
inclined  characteristic  1-3,  i.  e. 


lim 

i 

- 

oo  oo 


Thus  the  equilibrium  limit  of  Equation  (30)  after  (n  -  1)  iterations  is: 

Ap* 


lim 

i  .  _ 

T-TT 

00  00 


(66) 


Substituting  Equation  (66)  into  (56)  gives: 


Aq>-‘) 


Aq'  as 


T  U 
CO  00 


oo 


(67) 


From  Equations  (55)  and  (67)  we  see  that,  in  the  equilibrium  limit,  Aq*'  ' 
[jind  therefore  Ap'^^0  is  independent  of  the  number  of  iterations  for  n  S  Z. 
This  behavior  is  a  property  of  the  linearized  characteristics  and  cannot  be 
expected  to  be  valid  for  the  general  method  of  characteristics. 

In  Figure  2,  it  is  noted  that  increasing  c.  /c  (i.  e.,  increasing  the 

xoo  p 

energy  associated  with  the  relaxation  phenomena)  results  in  a  decrease  in 
Ap/Ap*.  Thus  the  effects  of  relaxation  on  the  wall  pressure  distribution 
increase  faster  with  temperature  than  the  error  introduced  by  the  finite 
difference  method  of  calculating  these  effects. 

If  the  equilibrium  characteristics,  Equation  (17),  were  used  to  cal¬ 
culate  pj  as  00  the  error  would  be  identically  zero  for  all  liet 

size.  Thus  there  is  some  large  value  of  i/r  U  at  which,  for  a  given 

"  oo  oo 

net  size,  it  would  be  better  to  calculate  the  flow  field  using  the  equilibrium 
characteristics.  In  Figure  4,  however,  it  is  seen  that  near  equilibrh'm  the 
error  is  strongly  dependent  on  the  location  of  the  initial  line  Th...  s 


due  to  the  dependence  of  the  error  on  the  gradients  of  the  thermodynamic 

properties  in  the  flow  field.  Thus  the  parameters  Ax/i  and  i/r  U  are 

00  00 

not  sufficient  to  establish  even  approximate  criteria  for  switching  over  to 

the  equilibrium  characteristics  as  i/r  U  becomes  large. 

00  00  ® 
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Figure  1.  Coordinate  System 
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